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Abstract 

In this paper we prove some characterizations of the matrix orthog- 
onal polynomials whose derivatives are also orthogonal, which gener- 
alize other known ones in the scalar case. In particular, we prove 
that the corresponding orthogonality matrix functional is character- 
ized by a Pearson-type equation with two matrix polynomials of de- 
gree not greater than 2 and 1. The proofs are given for a general 
sequence of matrix orthogonal polynomials, not necessarily associ- 
ated with an hcrmitian functional. However, we give several exam- 
ples of non-diagonalizable positive definite weight matrices satisfying 
a Pearson-type equation, which show that the previous results are non- 
trivial even in the positive definite case. 

A detailed analysis is made for the class of matrix functionals which 
satisfy a Pearson-type equation whose polynomial of degree not greater 
than 2 is scalar. We characterize the Pearson-type equations of this 
kind that yield a sequence of matrix orthogonal polynomials, and we 
prove that these matrix orthogonal polynomials satisfy a second order 
differential equation even in the non-hcrmitian case. Finally, we prove 
and improve a conjecture of Duran and Griinbaum concerning the triv- 
iality of this class in the positive definite case, while some examples 
show the non-triviality for hermitian functionals which are not positive 
definite. 
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1 Introduction 



The results published by Duran in jlUj can be considered the starting point 
for a general study of matrix orthogonal polynomials satisfying differential 
equations. After [TO], many other papers on the subject have appeared try- 
ing to find the similarities and main differences with respect to the classical 
and semi-classical scalar orthogonal polynomials (see jH [31 IHl El 1^ ITU 
I16j). In spit of these efforts, a complete Bochner-type classification of matrix 
orthogonal polynomials satisfying second order differential equations similar 
to the scalar case (see 012) is far from being obtained. 

However, the are many other differential properties that characterize the 
classical scalar orthogonal polynomials and that could lead to interesting 
matrix generalizations. These generalizations could clarify the structure of 
certain families of matrix orthogonal polynomials, being a source of proper- 
ties for such families, as in the scalar case. Eventually, the understanding of 
these other differential properties could shed light on the structure of some 
families of matrix orthogonal polynomials satisfying differential equations, 
helping to find classification theorems. 

It is well known that, apart from the second order differential equation, 
the classical scalar orthogonal polynomials {Pn) can be characterized by the 
orthogonality of their derivatives (Pn+i) (see P l8l ITTl [Tfll l2?1] ^ or, equiva- 
lently, by a linear relation between P„ and P^_^_l, Pl^, P^_i (see ^H])- Also, 
these properties are equivalent to a Pearson-type equation for the corre- 
sponding orthogonality functional (see [51 1191 I22j ). The main objective 
of this paper is to prove that the equivalence between these three properties 
hold in the matrix case too (see Theorem I3.14() . 

The proofs of the above equivalences are given for any sequence of ma- 
trix orthogonal polynomials, not necessarily related to an hermitian weight 
matrix. Consequently, the Pearson-type equation must involve a distribu- 
tional derivative. The distributional definition of the derivative not only 
permits to prove the results in a more general context, but unifies many 
different situations that otherwise would require a separate discussion. The 
reason is that the distributional Pearson-type equation takes care, not only 
of the first order differential equation for the weight, but of the necessary 
additional boundary conditions too (see Remark |2.9|) . So, the introduction 
of the distributional derivative becomes an advantage that permits to obtain 
more general results and, at the same time, in a simpler and more elegant 
way. 

Diagonalizable matrix orthogonal polynomials (we will be more precise 
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about this concept later) are nothing really different from scalar orthogonal 
polynomials. So, the relevance of the results proved in this paper depends 
on the existence of non-diagonalizable examples of matrix orthogonal poly- 
nomials whose derivatives are also orthogonal. Examples 2, 3 and 4 show 
that there are non-diagonalizable positive definite weight matrices whose 
orthogonal polynomials enjoy such a property. 
The weight matrix given in Example 2 



appeared previously in ^1] as an archetype of positive definite weight ma- 
trices whose orthogonal polynomials satisfy a second order differential equa- 
tion. Curiously, the authors declare in Jl], Section 7, Proposition 7.3, that 
the derivatives of these matrix orthogonal polynomials are no longer orthog- 
onal with respect to any weight matrix, arguing that a contradiction appears 
when supposing a three term recurrence relation for such derivatives. How- 
ever, if one makes the proposed computations in Proposition 7.3, no 
contradiction appears! Indeed, we will see that this weight matrix satisfies 
a Pearson-type equation that, according to Theorem I3.14( implies the or- 
thogonality of the derivatives of its orthogonal polynomials. Even more, we 
will find the positive definite weight matrix that gives the orthogonality of 
these derivatives. 

The purpose of Section 7, was to show that the equivalent charac- 
terizations of the classical scalar orthogonal polynomials do not necessarily 
hold for matrix orthogonal polynomials satisfying second order differential 
equations. It seems that the authors were not too lucky in the choice of the 
weight matrix since, if they had chosen the other example that they present, 
namely. 



they would have succeeded. The reason is that, as can be easily checked, this 
other weight does not satisfy the required Pearson-type equation and, then. 
Theorem 13. 141 implies that the derivatives of its orthogonal polynomials can 
not be orthogonal. 

A particular class of the family of matrix orthogonal polynomials with 
orthogonal derivatives permits a deeper analysis. This is the class corre- 
sponding to a Pearson-type equation involving a scalar polynomial a under 
the derivative. These matrix orthogonal polynomials can be classified anal- 
ogously to the classical scalar case, according to the roots of a: Hermite (no 
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roots), Laguerre (a simple root), Jacobi (two different roots) or Bessel-type 
(a double root). Moreover, a change of variable can reduce the different 
types to the canonical cases a{x) = l,x,l — x^^x^. 

For this special class we develop explicit formulas for the related ma- 
trix parameters, such as the norm of the monic orthogonal polynomials, the 
coefficients of the three term recurrence relation or the coefficients of the lin- 
ear relation between the polynomials and their derivatives. These formulas, 
although generalizations of the known ones in the classical scalar case, are 
more intricate due to the non-commutativity of the matrix product. How- 
ever, they are very useful since they allow to characterize the Pearson-type 
equations that have a quasi-definite solution. In other words, if a matrix 
functional satisfies this kind of Pearson-type equation, we have a criterion 
to know if it generates a sequence of orthogonal polynomials (see Theorem 
14. Notice that the importance of this result relies on the fact that we are 
dealing with general matrix functionals and not only with positive definite 
weight matrices, since the last ones always have an associated sequence of 
matrix orthogonal polynomials. 

We also prove that the matrix orthogonal polynomials of the above class 
satisfy a second order differential equation with polynomial coefficients (see 
Theorems 14 . 1^1 and 14 . 4j) . The result is again true no matter if the correspond- 
ing orthogonality matrix functional is hermitian or not. This is one of the 
novelties of this result, since the previous works on differential equations for 
matrix orthogonal polynomials always dealt with the hermitian case only. 
Indeed, if we believe a conjecture formulated by Duran and Griinbaum in 
jl3j , this discovering is only relevant for the functionals of the referred class 
that are not positive definite. This conjecture says that every positive defi- 
nite weight matrix in this class is diagonalizable. We present a proof of this 
conjecture (see Corollarv I4.11() . 

The above conjecture was supported on a partial proof given in _13_, that 
was incomplete due to the strong assumptions made there. First of all, it 
was supposed that the coefficients of the matrix polynomial appearing in the 
Pearson-type equation commute. Second, the proof was given separately for 
each of the canonical types of hermitian weight matrices that in the scalar 
case are positive definite: a{x) = — x"^. So, the case a{x) = is 

not considered, although the authors do not prove its incompatibility with a 
positive definite weight in the matrix case too. Finally, there is another less 
evident inconvenient. If a has a complex root, the required change of variable 
to arrive at a canonical situation destroys in general the hermiticity of the 
weight matrix. This means that, apart from the previous restrictions, the 
proof is only valid for the case of a with real roots. Our proof avoid all these 
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problems. Even more, we get a result that improves the one conjectured in 
jl3j (see Theorem I4.1U() . In spite of this result, the non-triviality of the 
class under consideration is ensured by the existence of non-diagonalizable 
matrix orthogonal polynomials in such a class, even in the hermitian case 
(see ini El and Example 5 of this paper) . 

The exposition of the above results will be structured in the following 
way along the paper. Section 2 introduces the notation, as well as some 
preliminary results and considerations that will of interest for the rest of 
the paper. In Section 3 we study the matrix orthogonal polynomials (P^) 
with respect to a functional satisfying a Pearson-type equation with two 
matrix polynomials of degree not greater than 2 and 1. We prove that such 
a Pearson-type equation is equivalent to the orthogonality of the derivatives 
(-P^+i) and, also, to a linear relation between P„ and P^_,_]^, P^, P^_^. Some 
two-dimensional non-diagonalizable examples of positive definite weight ma- 
trices whose orthogonal polynomials satisfy these properties are presented at 
the end of the section. Section 4 is devoted to the analysis of the special case 
in which the polynomial under the derivative in the Pearson-type equation is 
a scalar one. We obtain the characterization of the Pearson-type equations 
of this kind with quasi-definite solutions, the differential equation for the re- 
lated matrix orthogonal polynomials and the proof of the Duran-Griinbaum 
conjecture, finishing with some non-diagonalizable examples. Finally, in 
Section 5 we discuss the relation of the above results with other ones in the 
literature about second order differential equations for matrix orthogonal 
polynomials. 



2 The Basics 



We start with some notations and a summary of basic results that we will 
use in the rest of the paper. 

In what follows, will be the set of complex vectors of m components 
and C(™'™) the set of m X m complex matrices. We shall denote by P^*") 
the C^'"'™) -left-module 




and by means of p("^)' the C('"'"^)-right-module Hom(p("^), C^™'"")) . P^T^ 
will be the subset of matrix polynomials of p(™) with degree not greater 
than n. In the scalar case (m = 1) we will just write P^^-* = P and Pn ^ = Pn- 
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For all P e P^"*) and u G P^"")' the duality bracket is defined by (P, u) = 

u (P) and it verifies the usual bilinear properties. 

For G N and u G P*^"*-* the linear functional ux'^I G P*^'"-' is given by 

{P,ux''I) = {x''P,u), 

where / denotes the m x m identity matrix. A linear extension gives the 
right-product uQ G P^"")' for u G P^"")', Q G P^""), with Q{x) = Ylk=o<lkx'', 
Qk G C^"*'"*), in the following way: 

n 

{P,uQ)=^{x'P,u)qk. 

k=0 

Similarly, the left-product Qu G P*^"^-* is defined by 

{P,Qu) = {PQ,u). 

Every functional u G p(™) induces a matrix inner product in P^™-) given 
by {P,Q)u = {P,uQ*), where Q*{x) = YJk=oQk^^ and ql is the adjoint 
matrix of qk- This matrix inner product enjoys the standard sesquilinear 
properties. The orthogonality with respect to u means the orthogonality 
with respect to this inner product. 

The functional u* is defined by 

{P,u*Q) = {Q\uP*)\ 

and we will say that u is an hermitian functional if n = u*. In this case 
(P, uP* ) is hermitian for any P G p(™) . We will say that an hermitian 
functional u is positive definite if {P,uP*) is positive definite for every P G 
p("») with det P 7^ 0. In what follows we denote this condition by n > 0. In 
the same way, for a positive definite matrix A we will write A> 0. 

We denote by fik = {x'^I, u) the k-th. moment with respect to u e . 
Given a sequence (A*fe)fe>o C^™'"*), there exists a unique u G P("^)' such 
that {x'^I, u) = fik- 

If G P^™-* has moments (ma:)a;>o ' ^^y that u is quasi-definite (or 
non-singular) if det A„ 7^ for n > 0, where A„ is the Hankel-block matrix 



/no Ml • • • l^n \ 
Hi ■ ■ ■ /^n+l 

\ A'n /^n-t-1 • • • A*2n / 
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Notice that u is hermitian if and only if /i„ = ^* for n > 0, or, equivalently, 
A„ = a; for n > 0. 

The interest of the quasi-definite functionals rehes on the following result 

(see puniin])- 

Theorem 2.1. u € p("^) is quasi-definite if and only if there exists a se- 
quence {Pn)n>Q of left orthogonal matrix polynomials with respect to u, that 
is: 

(i) P„ gP^'"), degP„ = n. 

(ii) The leading coefficient of Pn is non- singular. 

(Hi) {x^Pn,u) = EnSnk, with En non-singular, for < k < n. 

Moreover, the sequence {Pn)n>o is unique up to non-singular left matrix 

factors and verifies a recurrence relation 

xPn (x) = OnPn+lix) + f3nPn{x) + 7„P„_l(x), 

where Pq G C*-™'™^ is non-singular, P_i = and an, Pn, 7n G C^™'™-*, with 
Oin, In non- singular. 

The last result of this theorem has a converse (Favard's Theorem): for 
any sequence {Pn)n>Q verifying the above recurrence relation there exists 
a unique (up to non-singular right matrix factors) quasi-definite functional 
u such that {Pn)n>o is its sequence of left orthogonal matrix polynomials 
(see jni I21j). Analogously we can define the right orthogonal matrix 
polynomials with respect to u, which are the adjoints of the left orthogonal 
polynomials associated with u* . In what follows we will consider only left 
orthogonal matrix polynomials, and we will call them just matrix orthogonal 
polynomials (MOP). 

Remark 2.2. Given a functional u G p(™)', we can normalize the correspond- 
ing MOP by choosing the only monic ones {Pn)n>o- In what follows we will 
assume this choice, so, a unique sequence of non-singular matrices {En)n>o, 
En = {x^Pn,u), is associated with any quasi-definite functional u. Also, 
and 7„ will denote the matrix coefficients of the related recurrence relation 

xPnix) = Pn+lix) + l3nPn{x) -F7„P„_i(x). 

Similarly, given a sequence MOP, we can normalize the corresponding 
functional u in different ways, for instance, by requiring (/, u) = I. How- 
ever, we will not fix the normalization for the moment because the most 
convenient one depends on the problem that we wont to study. 
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In the case of non quasi-definite functionals, the full sequence of MOP 
does not exist. Nevertheless, we have the following general result. 

Proposition 2.3. For every u G P^"*^ the following statements are equiva- 
lent: 

(i) Aq, . . . , A„ are non-singular. 

(a) There exists a finite segment (Pjfc)^=o ^/ iconic MOP with respect to u, 
that is: 

ra;PfcGP(™), dcgPfe = A:. 

(b) {x^Pk,u) = Ek6kj, with Ek is non-singular, for < j < k < n. 

Moreover, under the above conditions, the segment {Pk)k=o unique and 
there exists a unique monic polynomial Pn+i whit deg Pn+i = n -\- 1 such 
that {x^ Pn+i,u) = for < j < n. 

Proof. Suppose that Aq, . . . , A„ are non-singular. If Pk{x) = Z^f=o 7r|'^''a;*, 
g then, {x^Pk,u) = SLo^^^'Vi+i- Choosing Trf^ = /, the 

system Yli=o ''^f^^ l^i+j = 0, j = 0, . . . , A; — 1, can be represented as 

{t^^q \ T^f\ TT^^-l ) ^k-l = - (Mfc) • • • ) /^2fc-l ) , 

which has a unique solution for /c = 0,l,...,n+l. 

On the other hand, E^. is non singular for k = 0,1, ... ,n. In fact, we 
have {x^ Pk, u) = E^Skj, j = 0, . . . , k, k = 0, . . . ,n, and, so, 

W\4'\---,^^klvl)^k = i0, 0,...,0,Ek). 

If Ek is singular, there exists v G C"* \ {0} such that v^E^ = 0. Hence, 

( vT4'\ v^Trf\ . . . , v^irfl,, ) A, = ( 0, 0, . . . , 0, ) , 

and this result contradicts the non-singularity of A^ for /c = 0, . . . , n. 

For the converse, let us suppose that there exists a finite segment {Pk)k=Q 
of MOP with respect to u with E^ = {x^Pkju). It is easy to see that the 
conditions {x^Qk,u) = E^Skj, j = 0, . . . k, where Qk £ P^"^\ ensures that 

Qk = Pk^ k = 0, . . . ,n. Writing Qk{x) = Yli=o''^i'^^^^ above assertion 
means that, for = 0, . . . , n, the system 

( ttP, vrf ), . . . , 7rf_\, vrf ) ) A, = ( 0, 0, . . . , 0, ) 

has a unique solution and, hence, A^ is non-singular. □ 
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Concerning the partial hermiticity of a functional, we have the following 
immediate result. 

Proposition 2.4. Let u G p('")'. // {pk)l=Q is a basis o/P^T^ A„ = A; if 
and only if {(pkjUPj))]^ j=o hermitian. 

In particular, if u has a finite segment (-Pfc)^^Q of MOP, 

A„ = A; ^ {Pk,uP*) = Ek6kj, Ek = El 0<j,k<n. 

The second assertion of the above proposition says that, when Aq, . . . , A„ 
are non-singular, the condition A„ = A* means that the finite segments of 
left and right orthogonal matrix polynomials are each one the hermitian 
adjoint of the other one. 

Also, for the hermitian positive definite functionals on P^™^ we have the 
following characterization. 

Proposition 2.5. Let u € P(™) . // (pfc)fc=o i^ ^ basis o/Pi™\ the following 
statements are equivalent: 

(i) A„ > 0. 

(ii) {{Pk,up*))lj^,>0. 

(Hi) u has a finite segment {Pk)k=o of MOP such that {Pk,uPj) = EkSkj 
with Ej. > for < j,k < n. 

(iv) (P, nP*) > for any P G pI"^ such that det P ^ 0. 



Proof. We only prove (i) (iv), since the rest of equivalences are immedi- 
ate. For any matrix polynomial P{x) = X]i=o^«-^*' ^« ^ ^(m.,m)^ k < n, 



{P,uP*) = {Ai...Ak)A, 



A* 



So, {P, uP*) is hermitian if A„ is hermitian. If u G C™, 

v*{P,uP*)v = {vl ...vD^A ■■. y v^ = A*v. (1) 

Then, if t; 7^ 0, det P 7^ implies f i 7^ for some i. So, equality gives 
v*{P,uP*)v > if A„ > 0. 

For the converse, if (P, uP*) is hermitian for P G P^™'* with detP 7^ 0, 
l^2k = {x'^I, ux^I) = ^2k k < n. Besides, H2k-i = l^2k~i for A: < n too, 
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due to the identity {{x^ + ^)I,u{x^ + ^)I) = fi2k + m-2 + 2/Li2fc-i. 
Therefore A„ = A;. 

Suppose (P, uP*) >0 for any P G P^T^ with det P 0. Let [vq ■ . .Vk) , 
Vi G C", with Wfc / and k < n. We can always find Ai G C*^™'™) such 
that A*Vk = Vi, Ak = I. The polynomial P{x) = Yli=o ^^^^ '-'^ 
det P 7^ 0. So, relation gives 

{v*o ...vDAki ■■. \ >0, iivk^O, k<n. 

\Vk/ 

This proves by induction that A„ > 0. □ 

Remark 2.6. Notice that, if n is an hermitian and positive definite functional, 
then it is quasi-definite. So, there exits the corresponding sequence (Pn)n>o 
of MOP with En hermitian and positive definite. 

Similarly to the scalar case, the positive definite matrix functionals are 
those ones given by 

{P,u) = j P{x)dM{x), (2) 

where dM is a positive definite weight matrix on M, that is, a positive definite 
matrix of measures supported on the real line (M(S') is positive semidefinite 
for any Borel set 5" C M) with finite moments f x'^dM(x), n > 0, and 
such that / P{x) dM{x) P{x)* is non-singular if det P ^ (see jHI ) • This 
is, for instance, the case of an absolutely continuous matrix of measures 
dM{x) = W{x) dx with finite moments, W{x) being semidefinite positive 
for any x G M and non-singular for infinitely many points of the real line. 

In what follows we will identify any m x m matrix dM of measures 
on C with finite moments (not necessarily hermitian), and the functional 
u G P("")' defined by ©. Thus, we will write u = dM for such a functional. 

A specially interesting family of matrix functionals is given by the func- 
tionals which satisfy a differential equation of Pearson- type (see (HISl)- The 
definition of this family requires the introduction of the derivative operator 
in the space p("^) , which is the linear operator DiP^*") — > p(™) such that 

{P,Du) = -{P',u). 

The equality D{u<^) = (Du)^ + u<P' holds for ah u G P^*")' and ^> G P(™). 
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Definition 2.7. Let u G P^"*) . We say that u E V or, equivalently, is a 
P-functional, if there exist * G P^"*), with det $ 7^ 0, such that 

D (u^) = (Pearson-type equation) 

If dcg ^ < p and deg ^ < q, we say that u € Vp^q or is a Pp^g-functional. 
In both cases we also say that the corresponding sequence of MOP belongs 
to the family V or Vp q respectively. 

Remark 2.8. The condition det $ 7^ is imposed to avoid any triviality of 
the definition, ensuring that it involves all the components Uij : P^"*^ ^ C of 
u = ('Wij)fj=o- Notice that 

det $ = <^ ^v = for some v G C^lx] \ {0}. 

In fact, if = for some v G C^lx] \ {0}, then = (adj ^)^v = (det ^)v. 
To see the converse, remember that every $ G P^™) can be factorized as 
$ = P$Q, with d G P^"*) diagonal and P,Q e P^"") invertible, that is, 
det P, det Q G C \ {0}. Therefore, det$ = implies det 4> = and, since 
$ is diagonal, ^vq = for some vq G C"* \ {0}, which gives = with 
V = Q-^vo G C™[x] \ {0}. 

Remark 2.9. The distributional definition of the derivative operator D im- 
plies that, in general, the Pearson-type equation involves, not only a re- 
lation between standard derivatives, but a boundary condition too. Con- 
sider, for instance, a functional u = W{x) dx, x G P, with W an ana- 
lytic matrix function on a regular curve P of the complex plane. Then, 
Du = W'{x) dx + W{x){6{x — a) — d{x — b)) dx, where a and b are the initial 
and end points of P respectively. So, if the curve is open, together with 
the equality (W^)' = W"i>, we need the boundary condition (PP<5)(a) = 
{W^){b) = to ensure the Pearson-type equation D{u^) = u^. The case 
of a closed curve does not need an additional boundary condition since 
we suppose that W is analytic on P. Moreover, in this case, the Pearson- 
type equation holds even if (W^)' ^ but (W^)' - is analytic 
on the region enclosed by P, due to Cauchy's Theorem. The Pearson-type 
equation can be satisfied if W is only analytic on P \ {a, 6} but the limits 
{W^){a+) := limt^t„(W$)(7(i)), {W^){b-) := limt^tr{W^){l{t)) exist, 
where 7: [io,ii] — > P is a parametrization of P, a = ^{to), b = j{ti). Then, 

D{u^) = {W^)'{x) dx + (VP$)(a+) 6{x -a)dx- (VP$)(6") 5{x - b) dx, 
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so, we get the Pearson-type equation adding to (W^)' = the boundary 
conditions 

(W$)(a+) = {W^){b-) closed curve, 

(PF$)(a+) = {W^){b-) = open curve. 

The distributional derivative not only unifies all these cases, but allows 
to consider more general situations, such as functionals defined by matrix 
measures supported on an arbitrary subset of the complex plane. 

If u G P("^) is a P-functional with a Pearson-type equation D {u^) = u^, 
then, for every U G F^*") , 

D {um) = u ($0' + ^n) . (3) 

Therefore, the set 

M{u) = {$ G P^"") I D{u^) = u*, ^ e p(™)} 

is a right-ideal of p(™), but it is not necessarily principal, because the eu- 
clidean division algorithm is not valid in P^™). This is an obstacle to find 
a canonical representative of M.{u) that could lead to a classification of 
P-functionals similarly to the scalar case. 

Notice that V = \Jp,g>o'^p,q^ ^i^d 'Pp,q ^ "Pp',?' if P < p' and q < q'. The 

set 

Mp,q{u) = G Pj,'") I D{u^) = u^, ^ G fJ")} 

is not an ideal of P^™), but a C(™'™)-right-submodule of pj,"'^ Although it 
is finitely generated, it is not cyclic in general, what means again a problem 
for finding a canonical representative of Mp,q{u). 

Example 1. Let us consider u G P^^^' given by 

A direct computation shows that u is a 'P3^2-functional with 

M3,2iu) = spanj.(2,2) I (1 ~ ^^)-^' ^(1 ~ 
generated by two elements. Indeed, if 

$(x) = (l-x2)Ai+x(l-x2)('° J]a2, a, Gd2'2), 



2 m 
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then D{u^) = with 

We can get cyclic modules for u by going down in the net i'Pp,q)p,g>Oi 
but there are two different ways to do it. From the previous result we obtain 

• u e V2,2 with M2,2{u) = spanj.(2,2) |(1 — x^)/ 

• ti G Ps,! with 7W3,i(ti) = spanj-{2,2) |(1 — x^) 
In fact, 

_L). 

This splitting shows clearly the problem of classification of "P-functionals. 
Moreover, we can not go down more than this in the net (7^p,q)p,g>o since 

M2,i{u) = M2,2{u) n M3^i{u) = spanj;.{2,2) |(1 ~ 
Mi,2{u) = M3,o{u) = Mo,3{u) = {0}, 

and, hence, u Vp^q for p + q < 3. 

Notice that the above problems of classification happen even for quasi- 
definite functionals since our example was positive definite. However, if we 
restrict our attention to quasi-definite functionals, there is a singular situa- 
tion. As we will prove later (see Theorem K-i.4|) . if Aq, Ai, A2 are non-singular 
for some u € 7^2,1; then M2,i{u) is cyclic. This implies that we can associate 
with each sequence of MOP in the family 7^2, 1 a canonical representative: 
the unique (up to non-singular right matrix factors) generator of M.2,i{u), 
u being the related orthogonality matrix functional. 

A way to solve the problem of classification of P-functionals uses the fact 
that A4{u) always has a non-trivial scalar representative. In fact, choosing 
Q = adj $ in © gives ^0, = (det which yields the following character- 
ization (see W, 5_). 



(-2x l)}- 



: 
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Proposition 2.10. The functional u € p(™) belongs to the family V if and 
only if there exist a G P \ {0} and ^ € P*-™^ such that 

D{ual) = u^. 

Notice that the set 

M{u) = {a G P I D{ual) = u'^, ^' G p(™)} 

is a non-trivial bilateral ideal of P, which is, therefore, principal. So, there 
exists an a G P \ {0}, unique up non-trivial factors in C, that is generator 
of A4{u). This scalar generator can be used to classify the 7^-functionals. 

Definition 2.11. Let u G P(™)' be a P-functional and let a G P \ {0} be a 

generator of A4{u). The class of u is s = maxjdega — 2,deg^ — 1}, where 
qj g -g g^^j^ ^Yiat D{ual) = W^. 

The interesting "P-functionals are those ones that have a sequence of 
MOP, that is, the quasi-definite P-functionals. These are called semi-classical 
functionals (see I^IHl)- As in the scalar case, the semi-classical functionals 
can be characterized by several differential properties of the corresponding 
MOP. 

Theorem 2.12. Let u G P^"")' be quasi- definite and let {Pn)n>o 
associated sequence of MOP. Then, the following statements are equivalent: 

(i) n G P. 

(ii) There exist a G P \ {0} and 9^"^ G C^""'"") such that 

dcg a 

a{x)P'^^i{x) = Q^J^^ Pn+j{x) (structure relation) 

j=-s 

with s > maxjdega — 2,0} independent of n and G^j ^ for some 
n > s. 

(Hi) There exist a eF \ {0}, b e¥ and A^"^ G C^"''") such that 

a(x)f^'(x) + 6(x)P;(x) = ± Kf^Pr^^.ix) mero-d^fferent^al 

^ -TV/ equation) 

k=—r 

with r > maxjdeg a — 2, deg 6 — 1} independent of n. 
We use the convention P^ = for k < 0. 
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Proof. See HIS]. □ 

Remark 2.13. Let us suppose that a 'P-functional u G p(™') satisfies a 
Pearson-type equation D{ual) = u^, a G P \ {0}, ^ G P^™), and let 
s = maxjdega — 2,deg^' — 1}. Then, the proofs given in |S] show that 
the structure relation appearing in Theorem 12.121 (ii) is satisfied for the 
same polynomial a and integer s. However, contrary to the scalar case, the 
differo-differential equation given in Theorem 12.121 ( in ) can not be ensured 
for a = a., r = s., but for a = and r = max{2dega — 2, 2s + 2} = 
max{2 deg a — 2, 2 deg > s. 

In the scalar case, the classical orthogonal polynomials can be charac- 
terized by a Pearson- type equation D{ua) = uj3, a G P2 \ {0}, /? G Pi, 
for the corresponding orthogonality functional u. When trying to generalize 
the concept of classical orthogonal polynomials to the matrix case using a 
Pearson-type equation, the following two possibilities appear: 

• Zero class: u G P^™') belongs to the zero class if it is semi-classical with 
class s = 0, that is, u is quasi-definite and there exist a G P2 \ {0}, 
^ G P['"\ such that D{ual) = . 

• Family V2X- u G P^™)' is a P2,i-functional, or belongs to the family 
V2,u if there exist $ G P^"'', ^ G P^"\ with det ^> / 0, such that 

The MOP associated with zero class functionals or quasi-definite 7^2,1- 
functionals can be considered as matrix generalizations of the classical scalar 
orthogonal polynomials. Notice that a quasi-definite P2,i-functional is al- 
ways semi-classical, but its class can be greater than zero. In fact, excepting 
the scalar case, the family of quasi-definite 'P2,i-functionals is strictly greater 
than the zero class, as can be seen in Examples 2, 3 and 4. Both, the family 
V2,i and the zero class, are interesting sets of matrix functionals since the 
related MOP inherit some of the properties that characterize the classical 
orthogonal polynomials in the scalar case. This will be shown in the follow- 
ing sections, which are devoted to the study of the family 7^2,1 and the zero 
class. 

Before doing that, we will comment some other questions of impor- 
tance for matrix orthogonal polynomials. As we have pointed out, a central 
concept for matrix functionals is the diagonalizability or, more generally, 
the reducibility. We say that a functional u G P^"^) is diagonal or block- 
diagonal if its moment sequence {fJ-n)n>o enjoys such a property. We write 
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u = it^-^^ © • • • n^*^) if /x„ = iin^ © • • • © ijl^^ , where (iin^ )„>o are the moments 
of 

To simplify the analysis of a matrix functional u G ^ the usual 

strategy is to connect it with a diagonal or block-diagonal one u G P^"*) 
through a relation that permits to translate the information from u to u. 
For instance, if u = TuS, with T,S € £i"^^rn) non-singular, we say that u 
is equivalent to u. In particular, when S = T* we say that u is congruent 
to u, while ii S = T* = way say that u is unitarily similar to u. 
Notice the difference with the terminology used by other authors, we prefer 
to preserve the usual one in Linear Algebra to avoid unnecessary confusion. 
A matrix functional is diagonalizablc or reducible by equivalence if it is 
equivalent to a diagonal or block-diagonal one respectively. We define in 
a similar way the diagonalizability or reducibility by congruence and the 
unitary diagonalizability or reducibility. 

A change of variable t{x) = ax + b, a £ C \ {0}, b £ C, can be used to 
relate matrix functionals too. Given u G P^™) we define ut G P(™)' by 

{P,ut) = {Pot,u), 

so that, if M = dM, then ut = d{M o t~^). Notice that, with this definition, 

{Du)t = {Dut)t'. 

The kind of relation that we use depends on the properties that we need 
to preserve. For example, the equivalence transformation and the change of 
variable keep invariant the quasi-definite character, any family Vp^q as well as 
the class of a P-functional (in fact, the MOP and the corresponding Pearson- 
type equations are trivially related by these transformations). This means 
that, concerning these properties, the only non-trivial matrix functionals are 
those ones that are not reducible by equivalence or change of variable. In 
particular, if we are going to study a characteristic of a functional u that 
only depends on such properties, then we can always use the normalization 
(I, u) = I since we can work, for example, with the equivalent functional 
u = ujjLQ^. Also, this allows when studying zero class functionals to restrict 
our attention to the canonical choices a{x) = 1 — of the scalar 

polynomial in the Pear son- type equation, due to the freedom in the change 
of variables. 

However, if we are interested in a characteristic that depends on the 
hermiticity or positive definiteness of u (or, more generally, on the hermitic- 
ity or positive definiteness of some moments ^„ or Hankel matrices A„) we 
must use congruence transformations and changes of variable with real co- 
eflficients. This is the reason to avoid using the canonical forms of the scalar 
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polynomial a when studying hermitian zero class functionals, unless we are 
sure that a has real roots. Also, the normalization (/, u) = I can be used, 
while preserving any hermiticity property of u, whenever ^uq > since, then, 
we can use the congruent functional u = L~^ii(L~^)*, where //q = LL* is 
the Cholesky factorization of /xq- 

3 The family p2,i 

The aim of this section is to study the differential properties of the MOP 
associated with 7-'2,i-functionals. The main result is Theorem 13.141 which 
shows that some characterizations of the classical scalar orthogonal poly- 
nomials remain true for the matrix family V2,i- Along the way to prove 
Theorem 13. 141 we will obtain a chain of results which have their own inter- 
est. 

We will start fixing some notations that we will need in the rest of the 
section. Let u € P^™^ be a P2,i-functional, that is, D(w<I>) = n^', where 

<^{x) = ipQ + LfiX + (P2x'^, ^'(x) = Vo + i^lX, with (pi,1pj G £{m,m) 

det <I> 7^ 0. The above Pearson- type equation is equivalent to 

n(/i„_i(/?o + AtnV'l + = -(/^nV'O + /in+lV'l), « > 0, (4) 

where {iJ,k)k>o are the moments of u and = 0. We denote 

(JIq Jli ...Jin 


The moments of u and u are related by 

Jin = ^J'nVO + ^^n+Wl + ^^n+2^2■, n > 0. (5) 

One of the characterizations of the classical scalar orthogonal polyno- 
mials is that they are the only sequences of orthogonal polynomials whose 
derivatives are also sequences of orthogonal polynomials. The following 
proposition is the starting point to prove a similar result for the family P2,i- 

Proposition 3.1. Let u be a 'P2,i-functional such that Aq, Ai, . . . , A„ are 
non singular. Then, the corresponding finite segment {Pk)k^Q of monic MOP 
satisfies 

{x^Piu)=Q, j = 0,...,k-2, k = 2,...,n, 
(x'^-ip^.S) = -Efc(Vi + (A:-l)(/.2), k = l,...,n. 
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Proof. From the distributional equation D{u^) = we have 
or, equivalently, 

-j{x^~'Pk,u^) - {x^Piu<^) = {x^Pk,u^), 
which, for j = 0, . . . ,k — 1, gives the result. □ 

Corollary 3.2. Under the conditions of Provosition W^ (Pj^)^^]^ is a finite 
segment of MOP with respect to u if and only if the matrix + (A; — 1)922 
is non-singular for k = 1, . . . ,n. 

The above corollary shows the interest in finding conditions that ensure 
the non-singularity of the matrices ipi-\-k(p2, k = 0,1,2 ... . The next lemmas 
study the relation between the non-singularity of Aj, j = 0,1, ... ,p, and 
ipi + kip2, k = 0,1, . . . ,q, for small values of p and q. They also inform about 
the non-singularity of A^, k = 0,1, ... ,q, a result of interest since, in the 
scalar case, u is quasi-definite for any classical functional u. 

Lemma 3.3. Let u be a 1^2,1 -functional with Aq, Ai, A2 non-singular. Then, 
ipi and Aq are non-singular. 

Proof. If ipi is singular, there exists v G C™'\{0} such that ^iv = 0. Relation 
(jlj) for n = gives ^oV'o + /^iV'i =0- The non-singularity of fiQ = Aq implies 
ipQV = 0. So, from Q we have 

Hn-WoV + fJ-n^lV + fJ.n+imi' = 0, n > 1, 

and, hence. 











A2 1 


ipiv 


l-l 













Also, {ipov,ifiv,if2v) / (0,0,0) because det <^ / 0. Now, we can conclude 
the singularity of A2, with contradicts the hypothesis. So, ipi is non-singular. 

On the other hand, the calculation of i^i gives Ei = ^2 — /^iMo ^/^i; which, 
according to Proposition [231 is non-singular because Ai is non singular too. 
From (jSJ for n = we get 

Jlo = fJ-Qipo + fJ-lipi + /i2</'2, 
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and © for n = 0, 1 gives 



fiQlpO + IXilpi = 0, /io930 + ^19^1 + ^29^2 = -(/^i^o + M2^l)- 

Therefore, 

is non-singular. □ 

As a first consequence, we obtain the following announced result. 

Theorem 3.4. If u ^ V2,i and Ao,Ai,A2 are non-singular, the C*^™'"^^- 
right-module M2,iiu) is cyclic. 

Proof. Let us suppose that D{u^^^) = W^^'^ with ^>(*) e P^™\ ^(*) G pj""^ 
for i = 1,2, and assume that det<I>^^^ 7^ 0. We are going to prove that 
$(2) = ${i)A, A G C("'™). Let W(x) = with V-? e C^™'*"). 

Since u satisfies the hypothesis of Lemma l3.3| ■i/'^^^ is non-singular. Hence, 
A = $(i)(V'i^VVf ^ - ^^^^ satisfies 



From (01) for n = 0, ip^''^ = -/ig ViV'f\ therefore, = 0. If A{x) = 

Aq + Aix + ^2X^ A, G C(™'™), we get UnAo + + /in+2^2 = for 
n > 0, which implies 

Since A2 is non-singular, A = and, thus, <I>'^^) = ip^\ □ 

Now, we are going to consider 'P2,i-functionals satisfying the hypothesis 
of Lemma b. 31 Li such a case we can write tpi = I without loss of generality 
because the Pearson-type equation can be written as D{u^^Ij^^) = u^'Vf^- 

Lemma 3.5. Let u be a V2,i-functional with non-singular for k = 
0,1,2,3. Then, 

(i) ipi and tpi -\- ip2 are non- singular, 
(a) Aq and Ai are non-singular. 
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(in) u is a V2,i-functional, that is, D{u^) = u'$, with <I>(x) = X]^=o'^«-^*' 
^(^) = Sj=o^i'^"'' where (pi,'4'j € C^'"'™-' and det $ ^ 0. Moreover, 
<l, can 6e chosen such that (p2 = i^i^V2 o-nd ipi = + 2(f2)- 

Proof. We will assume without of loss of generality that ipi = I. 

(i) Let us suppose that I + ip2 is singular. There exists v S \ {0} such 
that ip2V = —V. Writing @ for n = 0, 1, 



m + noipo = 0, m^ipo + ipi)v + noifov = 0. 



Then, 



- IpoilpO + ipi)v + ifoV = 0. 
Consider @ again, but for n and n + 1: 

nun^iifo + Hn{ipo + rupi) + /U„+i(/ + nLp2) = 0, 
(n + l)/i„(/?o + Ain+i[V'o + (^^ + + p.n+2[I + {n + 1)992] 



(6) 



0. 



Multiplying the first equation on the right by "00 + Vi ^'^d subtracting the 
second one, gives 

n/i„_iv?o(V'o + V'l) + P-n [V'o (V'o + 931) - </?o + n (9^1 (V'o + 9^1) - m)] + 

+ [v92(V'o + V^i) - 9^1] - Ain+2 [/ + (n + 1)992] = 0. 

Then, taking into account ©, we get 

Ain-lV^O (V'O + ^l)v + Hn [ifl {ipO + Vl) " Vo] V+ 

+ Hn+l[^2{.'4'0 + ^l) - ^l]v - IJ'n+2V = 0, n > 1, (7) 



which implies 



A, 



[¥'i(V'o + ^i) - 'Polv 

— V 






Vo/ 



This contradicts the non-singularity of A3. 

(a) By Proposition 13.11 and Corollary 13. 2( {P[,P2} is a finite segment of 
MOP with respect to u. The result follows from Proposition 12.31 
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(in) The existence of matrix polynomials satisfying D(u^) = is 
ensured if 

+ = (8) 

Writing <l>(x) = ipo + Lpix + 'f>2x'^, ^(x) = + i^ix, (jHl) is equivalent to the 
system 








¥^0 






/ ^0 \ 






/ 


V'o 













^o(V'i - 2932) 





/ 















Vo 








1 + 2^2) 




\ ^2 ) 







(9) 



A solution of the last equation is ^1 = / + 2ip2, 922 = ¥'2- With this choice, 
converting the system into triangular form gives 

^ V'o ^1 \ ( \ ( \ 

1 (f2 ]\^^~]^\ 

^0 - 1pO(pi + i^o>f2 / \(Pl - i^Q / \-1po{VO-1p0^1+^l^2) ) 

From for n = 0, /xqV'o + /^i = 0, so, 

T := (^0 - V'o'/'i + ^1^2 = m + ^(TVi^^i + (^0 Vi)^'^2 = 

Since Ei = ^2 - fJ-il^o^ fJ'i, 

that, keeping in mind (@J for n = 1, can be expressed as 
T = -^(^'^(^iV'o + ^2 + -E'iy'2) = 

= -/^o ^(-^i/Ug Vl + Ai2 + -E'iy'2) = -fJ.Q^Ei{I + (/32)- 

That is, T is non-singular, what ensures that @ has a solution. 

Finally, we are going to prove that det <1> 7^ 0. From (jHl we can deduce 

- I.') = 

Since det <I> 7^ 0, det<I> = implies det(^' — <!>') = 0. However, taking into 
account that ipi = I + 2ip2 and (^2 = ^2 we get ^(x) — <&'(x) = V'o — + 
which has non-null determinant. □ 

Lemma 3.6. Let u be a V2,i-functional with non-singular for k = 
0,1,2,3,4. Then, 
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(i) ipi + j(p2 is non-singular for j = 0, 1, 2. 

(ii) Aj is non-singular for j = 0, 1, 2. 

Proof. We will assume without of loss of generality that V'l = I- 

(i) Taking into account Lemma 13.51 (Hi), the functional u satisfies D{u^) = 
, with ^2 = ¥^2 5 ^1 = + 2(/?2) where (pi, ipj have the same meaning as 

in the proof of the previous lemma. 

Let us suppose that / + 2(/P2_is singular. Then, there exists v G \ {0} 
such that cp2V = — ^v, that is, ipiv = 0. Since D{u^) = n^, we have 

niJln^l^O + Jj-n^l + V'n+W2) = -(fin'^Po + A'n+lV'l), « > 0, 

which, for n = 0, gives /IqV'o + /^iV'i = 0- Hence, ipQV = because, from 
Lemma 13.31 JIq = Aq is non-singular. So, 

(jln-l^O + Jj-n^l + Jln+1^2)v = 0, n > I. (10) 

According to ©, 

IJ-n-WO^OV + fInifl'fO + VO'fl)v + ^„+i((^2^0 + <^1^1 + V0'P2)v-^ 

+ /i„+2(<^2^1 + ¥'1^2)1' + lJ'n+3^2^2V = 0, n > 1, 

and from here we can deduce the singularity of A4, because ^P2^2V = ^l^v = 
^ 0. This contradicts the hypothesis. So, ipi is non-singular. 

(ii) From Corollarv l3.21 -P21 -^3} ^ finite segment of MOP with respect 
to u and, so. Proposition 12.31 ensures that A2 is non-singular. □ 

The previous lemmas can be generalized through an inductive process. 
This process will need the following result too. 

Lemma 3.7. Let u G p("")' and F G p{,"'\ with detF ^ 0. We denote 
u = uF and we suppose that there exist vo,vi, . . . ,Vq G C™, with Vk ^ for 
some k G {0, 1, . . . , q}, such that the moments (/in)n>o of the functional u 
satisfy 

^ fin+jVj =0, Vn > 0. 
i=o 

Then, there exist wi,W2,---,Wp-^-g G C", with Wk 7^ for some k G 
{0, . . . ,p + q}, such that the moments (/x„)„>o of the functional u satisfy 

p+q 

fJ-n+kWk = 0, Vn > 0. 

k=0 
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Proof. We will write F(x) = /q + /ix + • • • + fpxP with fi G C^™'™). Then, 
/i„ = X^f^Q fJ-n+ifi and the hypothesis of the lemma gives 

q q / p \ p+q p 

= ^ Jln+jVj = ^{^ f^n+j+ifi ]vj=^ fin+k ^ fiVk-i, 
j=Q j=0 \i=0 / k=0 i=0 

with the convention v-i = ■ ■ ■ = v-^p = 0. So, the vectors Wk = Yl^=o fi^k-i, 
k = 0, . . . ,p + q, satisfy the equality of the statement. It will be enough 
to prove that not all the vectors Wk are null. If all of them are zero, 
Yli=o fi'^k-i = iov k = 0, . . . ,p + q, and this implies 

p+q p 

k=0 i=0 

or, equivalently, 

<? / p \ 9 
= ^ ^ fix' Vj = F{x)J2 ^i^^ ' Vx G C. 

j=0 \i=0 / j=0 

Since det F 7^ 0, we obtain from Remark 12.81 that Yl'j=o ^i^"' ~ ^ 

X G C, which means that vj = for j = 0, . . . ,q, in contradiction with the 

hypothesis. □ 

Now we can reach the generalization of Lemmas 13.31 13.51 and 13.61 

Theorem 3.8. Let u be a V2,i-functional with non-singular for k = 
0,1, ... ,n, where n > 2. Then, Tpi + jip2 and Aj are non-singular for j = 
0,l,...,n-2. 

Proof. Due to Lemmas 13.31 13.51 and 13.61 the result is true for n = 2,3,4. We 
will assume the statement for an index n > 2, and we will prove that it is 
also true for n + 1. 

Assume that Aq, Ai, . . . , A„, A„_|_i are non-singular. Then, the hypoth- 
esis of induction implies that V'l + jV'2 and Aj are non-singular for j = 
0, 1, . . . , n — 2. We only must prove that tpi + {n — l)ip2 and A„_i are non- 
singular too. For this purpose we will introduce a set of 'P2,i-functionals 
= 0, 1, . . . , using the superscript {j) for the associated elements. 

Let us define n'-''-* = u, ^'^^^ = '^^^^ = ^ . Taking into account Lem- 
mas EH and given = nW$W(V'fV^ there exist G P^"'^ 
^(1) ^ pM^ satisfying D (tt(i)^>(i)) = ^(i)^'^!), with det$(^) / 0, = 
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Lf^'^ and ij)'^'^ = ijjf^ + 2ip^'^ non-singular. Moreover, from Proposition 13.11 
^(1) = + This implies that . . . , ^^^2 and, thus, 

Aq^\ . . . , A^"'^2 ai^6 non-singular. 

Following this procedure, we can construct inductively a set of 7^2,1- 
functionals u'^^\ j = 0, 1, — 1 {I = [§] ) , satisfying 

^(i+i) 

Aq-' A|j^22j non-singular. 

Let us suppose that n is even (n = 2/). Then, Aq ^\^2 
are non-singular. If -|- (n — 1)922 = + ¥^2 singular, the same 

arguments that lead to (4) in the proof of Lemma 13.51 give now 

3 

3=0 

Since u^^~^^ = uF, degF < 2/ — 2 = n — 2, we get from Lemma 13.71 

n+l 

fik+jWj = 0, some 7^ 0, k > 0. 

j=0 

This contradicts the non-singularity of A„+i, so, ipi + {n — l)(p2 must be 
non-singular. 

If, on the contrary, n is odd {n = 21 + 1), aJ"^\ aJ' A^'"^\ aJ 
are non-singular. Thus, analogously to (8) in the proof of Lemma 13.61 we 
find that, if ipi + {n — 1)992 = + 2(/?2 singular, 

4 

i=o 

Now, u^^~^^ = uF, degF < 2Z — 2 = n — 3, so. Lemma 13.71 gives again the 
same condition 

n+l 

fik+jWj = 0, some 7^ 0, ^ > 0, 

j=0 
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so, ipi + {n — l)ip2 is also non-singular in this case. 

Finally, the non-singularity of A„_i follows from Proposition 12 . 31 and the 
relation En-i = —^En{ipi + (n — 1)932) given in Proposition 13.11 □ 

The previous theorem and Corollary 13.21 have the following immediate 
consequences. 

Corollary 3.9. If u is a quasi-definite V2,i-functional, then ipi + n(p2 is 
non-singular for n = 0, 1, 2, . . .. 

Corollary 3.10. If u is a quasi-definite V2,i- functional, then u = is a 
quasi-definite V2,i- functional too. Moreover, if {Pn)n>o is the sequence of 
monic MOP with respect to u, then {^Pn)n>i ^"^ sequence of monic MOP 
with respect to u. 

Remark 3.11. The Pearson-type equation D{u^) = u^,<^ e P^"*\ ^ G P^™\ 
is equivalent to the recurrence nnn-i^Q+ tin{ipo+n^pi)+ fin+i{ipi+nip2) = 0, 
n > 0. Therefore, the non-singularity of the matrices ipi + nip2 for n > 
is a sufficient condition for the existence of a solution u of the Pearson-type 
equation. Indeed, this condition ensures that the solutions are determined 
by ^0 = (-^1 u) or, in other words, the solution is unique up to left matrix 
factors. Then, according to Corollary 13.91 if the Pearson-type equation 
has a quasi-definite solution, the quasi-definite solutions are exactly those 
solutions determined by a non-singular matrix fiQ. 



3.1 Characterization of the family 7^2,1 

In the scalar case, the classical orthogonal polynomials can be characterized 
alternatively by a Pearson- type equation (see 8j 19j ,20, jl2] ) ? the orthogonal- 
ity of the derivatives (see El I17| I19[ I2f)j ) or a linear relation between the 
polynomials P„ and P^+i, Pn-> ^n-i (see ^H])- The consequences of the pre- 
vious analysis provide an analogue of these equivalences for the matrix case, 
which constitute a characterization of the quasi-definite 'P2,i-functionals. In 
the proof of this characterization we will need the following results too. 

Lemma 3.12. Let u G P^™)' such that A n is non-singular. Then, 

uP = 0, PGP(;") P = 0. 

Proof Let P{x) = EiLo^^^*' ^' ^ C^™'™). Then, uP = is equivalent to 
UkAo H h /^fc+n^n — for A; > 0, which implies 
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and, thus, P = if A„ is non-singular. □ 

Proposition 3.13. Let u,v G ^uith u quasi-definite and (Pn) its corre- 

sponding sequence of monic MOP. Then, the following statements are equiv- 
alent: 

(t)v = uA, AeF^rl 

(a) (Pn) is quasi- orthogonal of order not greater than p with respect to v: 
{x''Pn,v) =0, k = 0, . . . ,n — p — 1. 

Proof. See ;5 . □ 

Here is the referred characterization of the quasi-definite P2,i-functionals. 

Theorem 3.14. Let u G P^*")' be quasi-definite and {Pn) its sequence of 
monic MOP. Then, the following assertions are equivalent: 

(i) u is a V2,i- functional. 

(a) (Pn) is a sequence of MOP with respect to a quasi-definite functional u. 
(Hi) There exist matrices an,bn € C^™'"^-* such that 

n+\ 

with jn — bn non-singular for n > 1. 
Moreover, n = n^>, ^> € P^"'\ det $ / and D{u^) = n*, ^' S P^™\ 
Besides, u is a quasi-definite V2,i- functional too. 

Proof. 

(ii) <^ (Hi) The sequence of matrix polynomials (P„) satisfies the recurrence 
relation, 

xPn = Pn+l + PnPn + InPn-l, 

SO, 

Pn = -XP'^ + P'n+l + PnPL + lnP'n-1- (H) 

If we assume (ii), {Pn) also satisfies a recurrence relation 

-XP^ = ^P'n+l + -~0n-lP'n + ^Tn-l^n-l (12) 

n n + \ n n — 1 



and, then, (fTT|) and (|T^ imply 

1 



Pn 



n + 



-P^+i + a„P; + 6„P^_i, (13) 
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where an = Pn-Pn-i and 5„ = 7n-;^7n-i- Notice that -fn-K = T^7n-i 
is non-singular. 

For the converse, from ((TT|) and ((T^. 

Now, we have a recurrence relation for (P^) with Pn-i = f^n — an and 
7„_i = ^^^^(7n — ^n)- Since 7^ — 6„ is non-singular, the Favard theorem 
assures the existence of a functional u £ P^™) such that (P^) is a sequence 
of MOP with respect to u. 

(a), (in) ^ (i) Assume the relation P„ = ;j^P^_|_i + o.nPn + bnP'n^i and 
the fact that (P^) is a sequence of MOP with respect to a certain functional 
u. Notice that this last hypothesis implies the non-singularity of En-i = 
^{x'^'^Pnju) for n > 1. Under the assumptions, 

{x''Pn,u) = 0, k = 0,...,n-3. 

So, (P„) is a quasi-orthogonal sequence with respect to u of order not greater 
than 2. Proposition 13. 13l savs that there exists ^ S Pg™^ such that u = u^. 
Setting w = D{u^), 

{x''Pn,w) = -k{x''-^Pn,u^) - {x^P'^,u^) =0, A; = 0, . . . , n - 2. 

Hence, (P„) is quasi-orthogonal with respect to w of order not greater than 
1 and, thus, there exists ^ G p(™') such that w = . 

It only remains to prove that det $ 7^ 0. For this purpose, notice that 
the equality 

{x^-^Pn,D{u<^)) = {x'^-^Pn.U^) 

gives 

-(n - 1)£'„9?2 - nEn-i = Enipi- 

Hence, ipi + {n — l)(/?2 is non-singular for n > 1. Suppose det<I> = 0. Then, 
according to R.emark 12.81 there exists v € C'"[x] \ {0} such that $z; = 0. 
Consider the matrix polynomial A G p(™) whose columns are all equal to v. 
Taking into account Lemma l3.12( the equality 

- ¥)A = (Du)^A = 

proves that {^! - (t')v = 0. So, "^v + ^-y' = and, if v{x) = H h Vnx'^, 

Vi E C", with Vn 7^ 0, we get {tpi + nip2)vn = 0, which is impossible. 
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(i) =^ (a) This impUcation is given by Corollary I3.1U1 □ 

Remark 3.15. Theorem 13. 141 ensures that any quasi-definite 'P2,i-functional 
u generates a sequence {u^"'^)n>o of quasi-definite ■P2,i-functionals, starting 
with u^^^ = u, and such that, for n > 0, 



u 



(n+l) ^ u^n)^{n)^ ^(n) ^ p^m)^ dct / 



Moreover, the A:-th derivatives {Prt'^)n>k form a sequence of MOP with 
respect to u^^^^ That is, as in the scalar case, if the first derivatives of a 
sequence of MOP are orthogonal, the higher order derivatives are orthogonal 
too. 

Remark 3.16. If u is not quasi-definite but Aq, . . . , A„ are non-singular, (ii) 
and (in) remain equivalent, but for the finite segment (Pfc)^^Q of monic 
MOP with respect to u. Besides, in this case, (i) also implies (ii) and (Hi), 
but only for the finite segment (-Pa;)^Zq , according to Theorem 13.81 

The following consequence of Theorem 13.141 will be of interest when 
studying the differential equation associated with the zero class MOP. 

Corollary 3.17. // a sequence (Pn) of monic MOP belongs to the family 
V2,i, then Pn±i G spaLn^m,^) {xP|^, P'^, Pn} . More precisely, 

where En = {x'^Pn,u), Pn{x) = x" + TrnX^~^ + • • • and Mn = Tpi + nip2- 

Proof. Using (|TT|) and (fTT?|) . we get by eliminating Pn+i and Pn-i respec- 
tively, 

r nPn = {x-pn + in+ l)an) PI, - (7n - (n + l)bn) 

1(1 - 6„7n') Pn=(^- bnln') P'n+1 + {bn%Hx - Pn) + ^n) 
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The matrix coefficients Pn^lm $mln,cim bn can be expressed in terms of En 
and TTn since 

/?n = 7I"n - 7I"n+l, In = EnE~^i, 

~ Pn—lj ~ Tn n—1 1 ' 

From fiere, it is just a matter calculation to get the result, using the fact 
that MfcM~^ = MfcMT^ = Mr^Mfe, where M„ = I + n(p2ipi^. □ 

3.2 Examples 

The purpose of the following examples is to show that non-diagonalizable 
matrix P2,i-functionals do exist, even in the positive definite case, and that 
the family V2,i is strictly bigger than the zero class (excepting the scalar 
case). Indeed, the presented examples are all positive definite and lie on 
the class s = 1. The matrix functionals of the examples have the structure 
w{x)R{x) dx, where is a classical scalar weight and 

p with positive leading coefficient, degg = l, 6 G C \ {0}. 

We will deal with a canonical form of these functionals, since any of them 
is congruent to one with the form 

(^^^''^"^J"!^''^ ""^^dx, pePmonic, aGC\{0}. 

This kind of functionals are never diagonalizable by congruence, neither by 
equivalence. This is a consequence of the fact that, as can be easily checked, 
any functional W{x) dx, with 



W 



Wn wi2 

W21 W22 



is non-diagonalizable by equivalence if {wii,wi2,W22} is linearly indepen- 
dent and {11^12,^^21} is linearly dependent. 

Example 2. Let u G P*^^)' given by 

u = e ^ ( ^ j dx, X eR, aeC\ {0}. 
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It is not a zero class functional, but its class is s = 1 due to the equality 

^ _ ( {\a?-2)x a \ 

^^-^^a(l-|a|V) -([a|2 + 2)xy' ■ 

Besides, it is a 7'2,1-functional with J^2,i{u) = spanj.(2, 2) {<!*}, where 

^{x) = -I 12 1 
The corresponding Pearson- type equation is D{u^) = n^', with 



-4x a 
2a -([a|2 + 2)x J ' 



Any right multiple of $ by a non-singular matrix factor can be chosen as a 
generator of A42,i (u), therefore, it will play a similar role in the Pearson-type 
equation for u. However, if we choose 



$(0) = $ 



1 
2 



the new functional u^^^ = u^^^^ is again a positive definite 'P2,i-functional 
of similar type. Indeed, 

(1) /|a|2 + 2-h2|a|22;2 2ax\ , 

^Uj = e ' ' ^_ ' ' ^ ] dx, X £R. 

\ 2ax 2 J 

This shows explicitly in the present example the general fact that any quasi- 
definite 7^2,1 -functional generates a sequence of 'P2,i-functionals, according 
to Theorem 13. 141 and Remark 13.151 

Example 3. The functional u € P*-^-* defined by 

x+|a|V aA^^ re G (0,00), a G C \ {0}, r > -1, 
ax 1 J 

lies again on the class s = 1 since 



, j r + 2 + {\a\'^ - l)x a 

^ ' * -a aPx^ r + I - {\aY' + l)x 



It is also a P2,i-functional, with M.2,i{u) = span£{2,2) {<!>} generated by 

^, , /(|a|2 + l)x 
V -a\a\ X X 
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The Pear son- type equation is D{u^) = -u*, where 



r + 2)a\a\'^x r + 1 - {\a\'^ + l)x J ' 
Notice that u^^^ = u$ is given by 

y ax 1 J 

so, it is a positive definite P2,i-functional of similar type. 
Example 4. The functional u G P^^^ given by 



u = x^e ^ 



is also in the class ,s = 1 since 



(r + |ap + 4)x — 

a(|ap + l)x2 (r - |a|^ + 2)x - 



and belongs to the family V2,i, with M.2,i{u) = span£.(2,2) {$} generated by 



X —a 

2 



(r + |a|^ + 2)xy ■ 
The Pearson- type equation is D{u^) = u'^, with 



(r + |ap + 3) - a; a 
-a(|ap + l)x (r + l)(r + 2) - (r + lap + 2)x ^ ' 



As in the previous cases, there is a choice of ^^'^^ € M2,i{u) that makes 
u^^^ = u<l>(^) a positive definite i-functional of similar type. The choice is 



and the new functional is then 



V (r + l)ax r + 2 J ' ^ ' ^ 
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4 The zero class 

The zero class is a specially simple subset of the family 7^2,1 • This simplicity 
allows for zero class functionals a deeper analysis than for general 7^2,1- 
functionals. According to the definition of the zero class we suppose in this 
section that u G p("^) is a quasi-definite functional that satisfies a Pearson- 
type equation 

D (ual) =u^, a G P2 \ {0}, ^ G P^^^ (14) 
We will use the notation a{x) = ao + Oiix + a2x'^, ai G C, and ^'(x) = 

The first aim of this section is to obtain explicit expressions for the 
elements associated with a zero class functional u in terms of the coefficients 
ai G C, il^j G C^™''™). This will lead to a characterization of the polynomials 

a G F2 \ {0}, ^' G pS""^ which can appear in the Pearson-type equation of 
a zero class functional. As a first restriction for a,'^, notice that Corollary 
13.91 implies that tpi + na2l must be non-singular for n > 0. 

Remember that (P„) denotes the sequence of monic MOP related to 
u, Pri{x) = x"I -|- TTriX^^"^ + ' ' ' and En = {x^Pn,u). As we have shown 
in the proof of Corollary I3.17| the coefficients of the recurrence xPn = 
Pn+i + PnPn + 1nPn-i and the Coefficients of the relation P„ = -|- 
anPn + bnPn-i cau be obtained from tt^ and En- So, we will just calculate 
TTn and En in terms of a and ^. 

From the Pearson-type equation for the functional u we obtain the rela- 
tion among the moments, that can be written in the following way 

nfJ.n~iao + llnNn + Mn+lM„ = 0, 71 > 0, (15) 

where Nn = V'o + nail., Mn = ipi + na2l- Taking n = and n = 1 in ()15|) 
we obtain 

W = -A^oV'o'0^^ H2 = fJ-o{'4'o'(Pi^'(po + aiipQtpi^ - ao)M{^. (16) 

Let us denote u = ual and (/i„)n>o its corresponding moment sequence. 
We know that 

/!„ = ao/"n + ai/"n+i + a2/^n+2, n > 0. 

This equality for n = 0, together with ((TB|). gives /io = fJ'00.{—ipoiPi^)'ipiM^^. 
Besides, a direct calculation shows that vri = — ^i/ig ^. So, 

^0 = Eoa{-i^o^^^)^iM^\ En = -^^„+iM„, 
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where = x" + + • • • and ^„ = n). 

Since u is a quasi-definite P2,i-functional, the same thing happens to u. 
Indeed, u is also zero class because D {ual) = n^', ^ = ^ + a'I. Notice that 
^(x) = V'o + "il^ix, where V'l = -^2 and ipo = Ni. 

The above results show that we can define a sequence {u^^^)j>o of zero 
class functionals by u^^^ = ua^ , and these functionals satisfy the Pearson- 
type equation 

D(n(^)a) = u'^^^¥^\ = '^ + ja'. 

Notice that Vo^^ = Nj, ivf,^'' = Nk+j, ip'f^ = M2j, Adjf^ = Mk+2j, where we 
denote with the superscript (j) the elements associated with the functional 
u^^\ Therefore, 

After an inductive process, 

vr. = vrf = nvri"-^) = uEf^Nr^.^M^^l^iEfY', 

En = = (-l)"n!4")M2-„i_2 • • • ^n-i = (-l)"n!4")M2„-iy-_\, 

where Vn = M„ • • • -^2n+i- Also, 

^J") = Eoa(-iVoMo-i)MoMf ^ • • • a(-iV„_iM2-J_2)M2„_2M2-^_„ 

and, so. 

En = {-irn\Eoai-NoM^^)MoM^^ ■ ■ ■ a{-Nn-iM^n^_^)M2n-2V-},. (17) 
If we define 11^ = E^^iTnEn, then 



(18) 



The above expressions give 7r„ and En in terms of a and ^ for a zero 
class functional u. When u satisfies the Pearson-type equation but it is not 
quasi-definite, the expressions for vTfc and Ek are valid for the finite segment 
{Pk)k=o of MOP with respect to u, whenever Aq, . . . , A„ and Mq, . . . , M2„-i 
are non-singular. This is because, then, the previous arguments remain valid 
for (n^-'^)"^Q and (-P^"'^)^"^, as follows from Corollarv 13.21 and Theorem 13.81 
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Moreover, if M2n, M2n+i are non-singular too, the formulas are also valid 
for the coefficients tTu+i, En+i of the extra polynomial Pn+i orthogonal to 
Pi"^\ given by Proposition 12.31 

With the achieved results we can get a characterization of the polyno- 
mials a, ^ related to the zero class. 

Theorem 4.1. The Pearson-type equation D{ual) = u'^ , a G P2 \ {0}, 
^ G p|^™\ has a quasi- definite solution u if and only if Mn and a{—NnM2n) 
are non- singular for n > 0, where Nn = ipo -\-naiI, Mn = ipi + na2l- Under 
these conditions, the solution of the Pearson-type equation is unique up to 
left matrix factors, and the quasi-definite solutions correspond to the non- 
singular choices of fiQ. 

Proof. If D (ual) = has a quasi-definite solution, the corresponding 
matrices En are non-singular for n > 0. Then, M„ and a{—NnM2n) are 
non-singular for n > 0, as can be seen from p7() . 

For the converse, from Remark 13. Ill if M„ is non-singular for n > 0, the 
solutions of the Pearson- type equation are determined by the choice of fiQ. 
Moreover, if, besides, a [—NnM2n) is non-singular for n > 0, the solution 
u is quasi-definite when /ig is non-singular. Indeed, proceeding by induction 
we can prove that there exist MOP with respect to u of any degree: 
■ There exists Pq = I, with Eq = fiQ non-singular. 

• Suppose that there exists a finite segment {Pk)k=o monic MOP with 
respect to u. By Proposition 12.31 there is a monic matrix polynomial Pn+i 
with degP„+i = n + 1, which is orthogonal to ¥^n^\ Since is non-singular 
for A; > 0, the expression of En+i = {x"^^^ Pn+i,u) is given by (|T7|l . Then, 
the non-singularity of a (— A'^fcM^^) for A; > implies that En+i is non- 
singular and, hence, {Pk)2t:0 also a finite segment of MOP with respect 
to u. □ 

Remark 4.2. From (|17() . we see that the non-singularity of for k < 2n — l 
and a{—NjM2j^) for j < n — 1, is equivalent to the existence of a finite 
segment {Pk)k=o of MOP with respect to any solution u of D (ual) = u\I' 
with fj,Q non-singular. 

As in the classical scalar case, every matrix functional in the zero class 
belongs, up to a change of variable, to one of the following types: 

• a{x) = 1 Hermite-type polynomials. 

• a(x) = x Laguerre-type polynomials. 
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• q(x) = 1 — Jacob i-type polynomials. 

• q(x) = Bessel-type polynomials. 

The characterization given by Theorem 14. II can be particularized for any of 
the above canonical types. For the Hermite-type polynomials, the existence 
of a sequence of MOP is equivalent to the non-singularity of ifji. In the 
Laguerre case, and ipQ + nl must be non-singular for n > 0. Jacobi-type 
polynomials exist if and only if tpi — nl and V'l i V'o ~ 2n/ are non-singular 
for n > 0, and, finally, the non-singularity of ipQ and V'l + f^I for > 
characterizes the existence of the corresponding Bessel-type polynomials. 
Notice that the conditions for the existence of Hermite, Laguerre, Jacobi 
and Bessel-type MOP are a natural generalization of the conditions in the 
scalar case. 

The non-singularity of the matrices M„ appeared previously in JTH' , as a 
condition for the Hermite, Laguerre and Jacobi-type polynomials to ensure 
that they are given by a Rodrigues formula. Our analysis proves that it is 
not necessary to impose this condition since it is automatically satisfied by 
any zero class functional. 

Theorem 14.11 has also important practical consequences for the study 
of MOP. When a matrix functional is given by a positive definite weight 
matrix on M, the corresponding MOP always exist. However, to decide if 
an arbitrary matrix of measures on M defines a quasi-definite functional can 
be a hard problem, even in the hermitian case. Theorem 14.11 solves this 
problem for any matrix functional satisfying a Pearson-type equation like 
(|14|) . Moreover, R,emark l4.2l gives a generalization that measures the length 
of the maximal finite segments of MOP associated with the functional when 
it is not quasi-definite. Some applications of this rule can be seen in Example 
5. The importance of the above result for the zero class will be clear later, 
since we will see that the only non-trivial matrix functionals in this class are 
not positive definite. 

4.1 Differential equation 

In this section we will prove that the MOP of the zero class satisfy a second 
order differential equation that generalize the known one in the scalar case. 
Notice that this is not ensured by Theorem 12.121 (in), since the right-hand 
side of the differo-differential equation given by this theorem could have 
more than one term, as follows from the comments in Remark |2^31 We will 
also obtain the structure relation of Theorem 12. 121 (ii). 
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In order to obtain the differential equation, starting from the study of 
the family 7^2, i) and keeping in mind Corollary 13.171 we can write for any 
sequence (-P„) of MOP in the zero class, 

p;±i = sWp„ + rWp^, (19) 

= (n + l)EnM^^^_,Mn^iE-\ 
Si"^ = -nEn-iM^^\^Mn-2E-\ 
' r^,+) = (n + l)EnM^^_^{a2E-^x - iM2„_2^-^7r„ + ^M2„_i£-i7r„+i), 
ri") = En-iM-^^_^Mn-2E-Hx + i7r„). 

On the other hand, Theorem l2.12l /^ii] and Remark l2.13l provide the structure 
relation 

aP^ = na2P„+i + r?nPn + 9nPn-l, Vn, On G C^'"'™). (20) 
Taking derivatives in the structure relation we obtain 

aP^ + a'P^ = na2Pn+l + VnPn + OnP'n-l 

and, using ((T^. we get 

aP^ + (a7 - r„)P^ - = 0, (21) 

r„ = na2rn'^^ + ^'nrn"^ + rjn, 
E„ = na2^n + ^n5^n , 

which is the differential equation for P„. 

We can calculate the coefficients of the above differential equation. First 
of all, notice that the coefficients r]n,On of the structure relation can be 
expressed in terms of 7r„ and En- A direct computation from the structure 
relation (|22|) gives 

r]n = nai + [(n - 1) vr^ - n7r„+i] 02, On = -EnMn-iE^\. 

Therefore, using 1)191) . (|21() and the above expressions, we find 

S„ = n^„M2~^_iM„_i[(n + l)a2 + Mn-2]En^ = nEnMn-iE-\ 
In the same way, writing r„(x) = Tn^x + F^f'*, Tn^ G ([;;;(™->™)^ -^^g gg^ 
r« = ^nM2-i_i [n(n + l)al - Mn-iMn-2] = -^„M_2^-\ 
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4°) = noi - -EnM-]^^ [n(n + l)a2M2n-2+ 

+ M„_iM„_2 - n{n - l)a2M2„-i] = 

= nai EnM2n-2E~^T^n = nai - -EnM2n-2^nE~^ , 

n n 
where Tin is given in (|18() . From (|18() and the above result we finaUy obtain 

a'{x)I - Tn{x) = En^iE-^x + EnVn-i'4>oV-\E;^\ 

Summarizing, we can enunciate the following result. 

Theorem 4.3. Let u he a zero class functional with Pearson-type equation 
D{ua) = n^-, a G P2 \ {0}, ^ G 

(i) If {Pn) is the unique sequence of monic MOP with respect to u, 

aPl[ + EnVn-i-^V-\E-^Pl, - nEnMn-iE-^Pn = 0, 

where Mn = 4^1 + 7102/ and Vn = MnMn+i ■ ■ ■ M2n+i- 

(ii) If {Qn) is the unique sequence of MOP whit respect to u such that Qn 
has a leading coefficient Kn = {EnVn-i)~^ , 

aQl + ^Q'^ - nMn-iQn = 0. 

The differential equation satisfied by the MOP of the zero class charac- 
terizes such MOP, as the next result shows. 

Theorem 4.4. Let u he a zero class functional with Pearson-type equation 
D (ual) = n^, a G P2 \ {0}, ^' G Pj^™\ Then, the differential equation 

ay" + W - nMn-iy = 

has a unique (up to right matrix factors) matrix polynomic solution y G p(™). 
This solution is the only n-th MOP Qn with respect to u which has a leading 
coefficient k„ = (E'nK-i)"^ • 

Proof. Trying y = Y2k>o ^kX^ as a solution of the differential equation, we 
obtain the recurrence for the coefficients 

(n - k)Mk+n~iCk = (k + l) [NkCk+i + (A; + 2)000^+2] • 

Since is non-singular for n > 0, for every k n, c^+i = Ck+2 = implies 
Cfc = 0. Hence, any non-trivial polynomic solution must have degree n, and 
such a solution is determined by c„. 

If Cfc = for k > n and Cn = Kn, there exists a unique solution that must 
be Qn- If, on the contrary, = for k > n but c„ is arbitrary, the solution 
is QnLn, where L„ = K~^Cn. □ 
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4.2 The hermit ian case 

Among all the zero class functionals, the hermitian ones have remarkable 
features that deserve to be emphasized. Maybe one of the most important 
has to do with the diagonalizatibility. 

The main purpose of this section is to prove a conjecture of Duran and 
Griinbaum (see ^Sl)- any positive definite zero class functional is diagonal- 
izable by congruence. In fact, we will prove a more general result, since we 
will get the unitary diagonalizability and, at the same time, under much 
weaker conditions for the matrix functional. The key result to prove the 
referred conjecture is the following one. 

Proposition 4.5. Let u G P^™)' be a solution of D (ual) = u^, a G P2\{0}, 
^ G p['"\ // fj.n-2, ■ ■ ■ , are hermitian, 

with Aq = ^(aoai), Ai = 29(aoa2)i A2 = 9(0102). 
Proof. From the hypothesis. 

Let us calculate {^*x"',u'^). 

(^*x",u^') = (^*x",L'(ua)) = -n(^*x"-\uo) - Vi*(x",uo) = 

= -n(Qx"-\n^)* - ((ax"-\n^) - (ax"-\ n)V'o)* = 
= -{n + l){ax''-^ , D{ua))* + ipoix""-^ ,ua) = 

= -(n + l)(ax''-\D(ua))* - -ijUx'',u^). 

n 

Using the above results we get 

{n + l){{ax"-~'^ , D{ua)) - (qx"'\ I)(ua))*) = i(-0oMn+i'0i - 

n 

which, together with the equality 

(ax"^S D{ua)) = -{n - l)(|a|2x"-2, u) - {a'ax'^-\u), 

gives 

ipQfin+iipi - V'l^n+iV'o = "-("- + l)((aa' - a'a)x""\ n) = 
= i2n{n + 1) [Q{aoai)fin-i + 29(aoa2)/in + 9(aia2)Ain+i] • 

□ 
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Using the standard notation [A, B] = AB — BA for the commutator of 
two square matrices A, B, we get the foUowing immediate consequence of 
Proposition 14.51 

Corollary 4.6. Under the conditions of Proposition if fiQ = I and 
is hermitian too, 

with the coefficients Aq^Ai, A2 as in Proposition \4-5\ 

The commutativity of a set of hermitian matrices is equivalent to state 
that they are simultaneously unitarily diagonalizable. Therefore, Corollary 
14.61 relates the possibility of diagonalizing simultaneously and /ii, to 
the requirement for a to have real coefficients. The next theorem gives 
conditions which ensure that a must be a real polynomial. 

Remember that, if ^0 > for a matrix functional, we can normalize it 
by /^o = without loosing any hermiticity property of the functional. So, in 
what follows, we will use freely this normalization when it is possible. 

Theorem 4.7. Let u G P^™)' he a solution of D (ual) = u^, a G P2 \ {0}, 
/u„ = /u* for n < 5, then a is a real polynomial (up non-trivial 
factors ) under any of the foUowings conditions: 

(i) [^2i /Ui] = 0, Aq > and Ai, . . . , A5 non-singular. 

(ii) A2 > 0. 

Proof. Without loss of generality, we can suppose //q = I- Let Aq,Ai, A2 be 
the coefficients given in Proposition 14.51 

(i) [Ei,fii\ = since Ei = ^2 — Mi- Then, from ((T7|) for n = 1, we obtain 
[V'i,A*i] = 0, which implies [V^Oj/^i] = because V'o = — /^iV'i- Using (fT5|) 
and the fact that Mn is non-singular for n < 3, due to Theorem 13.81 we get 

= for n < 4. Then, from Corollarv 14. 6| 

which implies = 0, Vi. 

(ii) Corollarv 14.61 for n = 1, 2, 3 gives 
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Therefore, 




(6^0 [1^2,1^1] + 2^1 [//3, 1^1] + M [M4, Ati]) V'l- 



Notice that, if P{x) = (Aq + + A2X^)1 
{P,uP*) = (^0 ^1 A 




Let us suppose P ^ 0. Since A2 > 0, Proposition 12.51 impUes that 
{P,uP*) > 0. From Lemma we know that Vi is non-singular, so, the 
matrix {■ip^^)*{P, uP*)ip^^ must be positive definite too. On the other hand, 
tr[^„,/ii] = and, thus, tr ((V'r^)*(P, MP*)V'r^) = 0- Hence, {P,uP*) can 
not be positive definite. This means that P = and Ai = 0, Vi. □ 

Corollary 4.8. For any positive definite zero class functional, the scalar 
polynomial of the Pearson-type equation is real up to non-trivial factors. 

The following result says that a zero class functional with a real scalar 
polynomial in the Pearson-type equation does not need to many conditions 
to be unitarily diagonalizable. 

Theorem 4.9. Let u be a zero class functional with jjLn = fi^ /''^ n < 3 and 
Aq > 0. Then, if the scalar polynomial a of the Pearson-type equation is 
real up to factors, u is unitarily diagonalizable. 

Under the above conditions, i/^4,//5 are hermitian too, then, a is real 
up to factors if and only if u is unitarily diagonalizable. 

Proof. Suppose without loss of generality that //q = I- If ^« = 0, Vz, 
Corollarv 14.61 for n = 1 gives Vi [/^2,A'i]V'i = 0. Since V'l is non-singular, 
[^2l^J'l\ = 0, so, there exists T € C^™'™^ unitary such that TfinT* is diago- 
nal for n = 1, 2. Then, TEiT* is diagonal because Ei = ^2 — fJ-i- From (|17|) 
for n = 1 we find that TipiT* is diagonal too. Hence, Ti^qT* is also diagonal 
due to the identity ipQ = —fiiipi. Using (|15j) and the non-singularity of M„ 
for n > one finds that TfinT* is diagonal for n > 0. 

The converse when ^4 , ^5 are hermitian follows from Theorem l4.7l (i). □ 

Joining Theorem 14.71 and 14.91 we achieve the following result that goes 
even further than the conjecture of Duran and Griinbaum. 
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Theorem 4.10. Let u be a zero class functional with fin = f^n f^''" n < 5. 
Then, u is unitarily diagonalizable under any of the followings conditions: 

(i) Aq > and [/i2,/Ui] = 0. 

(ii) As > 0. 

Notice that some of the conditions in Theorems 14.71 14.91 and 14.101 can 
be weakened. For example, in Theorem 14.71 (i), it is possible to substitute 
the condition Ai,...,A5 non-singular by A2 non-singular and [/i3,;Ui] = 
[M4, W] = 0. 

Corollary 4.11. (Durdn-Griinbaum conjecture) Any positive definite zero 
class functional is unitarily diagonalizable. 

The above result does not mean that the hermitian zero class is trivial, 
since there exist non-diagonalizable zero class MOP with respect to her- 
mitian functionals which are not positive definite (see Example 5). What 
is trivial is the positive definite subclass of the zero class (actually, a big- 
ger subclass, according to Theorem 14. 10(1 . Hence, positive definite Hermite, 
Laguerre and Jacobi-type MOP are unitarily diagonalizable. Concerning 
the Bessel case we can say even something more: similarly to the scalar 
situation, positive definite Bessel- type MOP do not exist, as the following 
proposition asserts. 

Proposition 4.12. Any zero class functional whose Pearson-type equation 
has a scalar polynomial with a double root, is not positive definite. 

Proof. Assume that u is a positive definite zero class functional whose cor- 
responding Pearson- type equation has a scalar polynomial a{x) = (x — a)^, 
a G C. From Corollarv 14.81 a G M. Also, Corollary 14.111 implies that there 
exists T G C^™'™"^ unitary such that TfinT* is diagonal for n > 0. There- 
fore, TEiT* is also diagonal and, using we get that TipiT* and TtpoT* 
are diagonal too. So, if we define the change of variable t{x) = x — a, the 
diagonal hermitian matrix functional ut = TutT* satisfies the Pearson-type 
equation D{utt^I) = utT^it + a)T* . Hence, ut = u't'^ • • • © u^^\ where 
u^f' are scalar Bessel functionals. Since a scalar Bessel functional can not be 
positive, the functional u is not positive definite, in contradiction with the 
hypothesis. □ 

4.3 Examples 

An example of non-diagonalizable hermitian zero class functional was pre- 
sented in [S]. JSI generalizes this example and provides several non-trivial 
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families of hermitian matrix functionals that satisfy a Pearson-type equa- 
tion like (|14() . In this section we will use these examples, including some 
non-hermitian generalizations, and we will prove that the corresponding 
zero class MOP do exist as an application of Theorem 14. II Notice that ^H] 
does not answer this question since the analysis of the non-positive definite 
weights dM{x) given there was under the assumption that P{x) dM{x)P{x) 
is non-singular for any matrix polynomial P with non-singular leading coef- 
ficient, something that was not proved in the concrete examples. 

The non-diagonalizability of the functionals given in the following exam- 
ples is ensured because they have the structure u = W{x) dx, where 

^r _ I Wu Wl2 



W21 

with {tL'ii,t«i2} linearly independent and {1^12, ^^'21} linearly dependent. 
These conditions imply that the functional u is not diagonalizable by con- 
gruence or, even, by equivalence. 

Example 5. Let us consider a functional u G P^^^ given hy u = w{x)R{x) dx, 
where w is a positive classical scalar weight with Pearson equation (wa)' = 
wfi and 



R{x) 



^ + geP,\{0}, a,6G C\{0}, CGC. 



Notice that u is hermitian when b = a, c G M and g is a real polynomial. 

This kind of functionals always satisfy the boundary conditions which 
ensure that D[ual) = {ual)' (see Remark 12. 9|) . In fact, writing them in the 
canonical representations, they have the form 

e-^(^ + ^^^^+"2^' lyx, XGM, 
^r^-.^c + c,x + c,\og{x) a^y^^ ^e(0,oo), 

(i+xr(i-xr(^+'^^^°s(i+x)+c2iog(i-x) a^y^^ xe{-i,i), 

in the Hermite, Laguerre and Jacobi case respectively. In the above expres- 
sions ci , C2 G C do not vanish simultaneously and r,s > — 1 . 
The functional u satisfies the Pearson-type equation 



D{ual) = W^, * 



P 
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Therefore, if q{x) = qo + qix and (5{x) = /3o + l^ix, 

Pi + na2 



21 

a 



a{-NnM^^) = a 



Pi + na2 



Po + nai \ f I 



Pi + 2na2 y V * 1 



Notice that, due to Theorem l4.11 Pi+na2 and a(— ^°_^2na^2 ) ™ust be different 
from zero for n > 0. Hence, M„ and a{—NnM2^) are non-singular for n > 0. 
Also, /Uo is non-singular since 



Mo = t'o(^^ o) ' ^0 = j^w{x)dx. 



So, according to Theorem 14. H we conclude that the functional u defines a 
sequence of zero class MOP. 

The above two-dimensional examples are only particular cases of the 
m-dimensional zero class functionals belonging to the equivalence classes 
defined by 

gAx^-Bx^^^^ X G M, K(A) > VA G spec(B), 



x\-^^dx xe(Ooc) I > -1 VA G spec(^), 
ax, xfc^u,ooj, <^ s)^(^) > VA G spec(B), 



{l + x)^{l-x)^dx, xG(-1,1), K(A) > -1 VA G spec(^),spec(5), 

where A,Be commute and spec(j4) means the spectrum of the 

matrix A. The restrictions for the spectra ensure the integr ability for any 
matrix polynomial and, together with the commutativity of A and B, lead to 
a Pearson-type equation of Hermite, Laguerre and Jacobi-type respectively, 
according to Remark 12.91 The conditions for the spectra also ensure the 
existence of MOP whenever is non-singular, as follows from Theorem l4.11 
For some choices of A and B it is possible to get an equivalent hermitian 
functional. This is the case of the initial examples, as 13 points out. 

These examples do not cover the zero class functionals of Bessel-type. 
Such examples can be found starting from a scalar Bessel weight. For in- 
stance, w{x) = x^e^/^, with r = —1,0,1,2,... , is a Bessel weight on the 
unit circle T := {x G C | |x| = 1} with Pearson equation (wa)' = wP, 
a(x) = x^, P{x) = {r + 2)x — 1. The matrix function W = wR satisfies 
the equation (Wa)' = W^, where R and ^I' have the same meaning as 
previously. However, 



W{x) = x^e 



r l/x 



C+ — +C2 log(x) 
X 

b 
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is not analytic on T if C2 ^ 0. If, for instance, we choose a logarithm with 
the discontinuity at the non-negative real axis, the matrix functional u = 
W{x) dx, X G T, verifies (see Remark 12. 9|) 

D{ual) = {Wa)'{x) dx — i27rec2 ^ ^ ^ ^ 5{x — 1) dx, 

so, it satisfies the Pearson- type equation D{ual) = when C2 = 0. 

Similarly to the initial examples, this new one is equivalent to a particular 
two-dimensional case of the general m-dimensional zero class functionals 
with the form x^e^^^dx, x G T, where r = —1,0, 1,2, .. . and B G C^'"''") 
is non-singular. Analogously to the scalar case, these functionals satisfy 
a Pearson-type equation of Bessel-type since the restriction on r gives the 
analyticity on T for x^e^^^. As in the previous examples, the conditions 
for r and B ensure the existence of the corresponding MOP when hq is 
non-singular, due to Theorem 14. 11 

Concerning the restriction on r it is known that, for the Bessel scalar 
case, it can be weakened to r ^ —2, —3, ... by introducing the alternative 
weight on T 

r(r + 2) 1 

This weight satisfies the equation (woa)' = wq(3 + r + 1, a{x) = x^, (3{x) = 
(r-f 2)x— 1. So, according to R,emark [2.9| the scalar functional uq = wo{x) dx, 
x € T, verifies the Pearson-type equation D{uQa) = tio/3. 

Notice that = t \\ where, in general, we denote 



r(r+2+fc) ^f+2y 



(I iik = 0, 

\A{A + I)---{A + {k-l)I) if ken, 



for any square matrix A. lfA,B£ C^™'"") and spec(y4)n{0, -1, -2,...} 
we can consider the matrix function 



oo ^ 



k=0 

which is analytical on C\{0}. If, besides, A and B commute, then (Wa)' = 
W"^ + A — I, a{x) = x^, ^'(x) = Ax — B. Hence, the matrix functional 
u = W{x) dx, X G T, satisfies the Pearson-type equation D{ual) = 
analogously to the scalar case. Therefore, Theorem 14.11 states that there 
exist Bessel-type MOP associated with u when B and /xq are non-singular. 
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5 Other differential equations 



Among the results proved by Duran in ^U], we remark in this section one 
concerning the existence of differential equations for MOP with respect to 
hermitian functionals u € P^™) satisfying a Pearson-type equation 

The referred result states that such a Pearson-type equation, together with 
the hermiticity of u^, is equivalent to state that the corresponding MOP 
(Pn) satisfy a second order differential equation 

P^'^>* + P^-f * + AnPn = 0, (22) 

with An G C^'"''"'* such that A„(P„,P„)„ is hermitian (actually, the re- 
sult is proved in jlOj for matrix orthonormal polynomials with respect to 
positive definite matrix functionals, but the generalization to the quasi- 
definite hermitian case is immediate). If, as in the rest of paper, we suppose 
that the MOP are monic, the condition for A^ becomes A^En — Pn^n- 



Also, equaling the coefficients of the highest powers of x in (|22j) we get 
An = -n{n — 1)^1 — rnp2 = -nM*_^. 

All the examples of 'P2,i-functionals u € P^'^^ presented in Section 3 
were hermitian and positive definite and, for all of them, we found a matrix 
polynomial <I> G ■M.2,i{u) with det<I> ^ such that is also hermitian 
and positive definite (in Examples 2 and 4 such a matrix polynomial was 
denoted 'I'^'^^ we omit now the superscript for convenience). Therefore, the 
corresponding MOP (P„) must satisfy a second order differential equation 
like 

For instance, in the case of the functional given in Example 2 



1 + laPx^ ax 



ax 



'^^dx, xGM, aGC\{0}. 



we find 



\(\(^?' + '^ -a\a\^x\ . f -4x 2a \ 

^"^^^V 2 )^^^^^'\2a -2(|a|2 + 2)xj + 

+ ^(o 2(jap + 2))^"(^) = 0- 

This functional was previously studied in |T4|, where it was proved that the 
corresponding MOP satisfy another second order differential equation lin- 
early independent with respect to this one. The fact that, contrary to the 
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scalar case, the MOP can satisfy linearly independent second order differ- 
ential equations was recently discovered (see [71 ITB]). 
As for the functional 



u = x^e ^ 



('^^Ix^''^ T)"^^' aeC\{0}, r>-l, 



given in Example 3, we get 



^ P'^^W (^ + 2)(|a|2 + l)-x -{r + 2)a\a\^x \ 
+ _ r+l-{\a\^ + l)x)^ 

(\ (r+ l + n)ala|2\ ^ , . ^ 
+ nL P„(x)=0. 



Finally, Example 4 deals with the functional 

ax 

whose MOP must satisfy the differential equation 



u = X e 



IX \ 

^ jdx, xe{0,oo), aGC\{0}, r > -1, 



P"(r) / + 1 + 



" Y a 



x] — (r + l)a([a|2 + l)x 

(r + l)(r + 2) - (r + jap + 2)x 



Let us restrict our attention now to the zero class MOP, that is, those 
whose corresponding functional u G P^™) satisfies a Pearson-type equation 



D{ual) = li^, a e P2 \ {0}, ^ G 



If It is hermitian, the hermiticity of ual is equivalent to saying that a is a 
real polynomial. Hence, if u is hermitian and a is real, the MOP (Pn) with 
respect to u satisfy the second order differential equation 

qP^' + P^^* - nM*_^Pn = 0. 

This differential equation is similar, but not equal to the one given in The- 
orem 14.31 However, when /ig > this difference disappears since, then, 
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Theorem 14.91 implies that u is unitarily diagonahzable. That is, there exists 
T G C^™'*"-* unitary such that u = TuT* is diagonal hermitian, so, the cor- 
responding monic MOP (P„) must be diagonal with real polynomials in the 
diagonal. Following similar arguments to those given in the proofs of the 
theorems in Section 4, we find that ^ = T^T* is also diagonal. Moreover, 
D{ual) = u^, hence, ^ is real. Therefore, both differential equations are 
the same for and, thus, also for (P„) since P„ = TP„T*. 

Returning to the family V2,i, the two-dimensional examples that we have 
found suggest that, for a big subclass of hermitian P2,i-functionals, the re- 
lated MOP satisfy a second order differential equation like (f^ . Equiva- 
lently, it seems that for many hermitian 7'2,1-functionals u E P^*") it is pos- 
sible to find a generator $ of the module A^2,i {u) such that is hermitian 
too. In particular, the referred examples seem to indicate that if u is positive 
definite, then is also positive definite for some generator $ of A42,i(w). 
The characterization of the subclasses of hermitian P2,i-functionals which 
are invariant under the operation u (for some choice of the generator 

$ of A^2,i('u)) remains as an open problem. This is an important ques- 
tion, not only for the study of differential equations for MOP, but also for 
the development of a general and systematic method to obtain modified 
Rodrigues' formulas for 'P2,i-functionals (see ^Ij for some examples of this 
kind of Rodrigues' formulas), as it will be shown in a future paper. 
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